Cellular Dynamical Mean Field Theory of the Periodic Anderson Model 
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We develop a cluster dynamical mean field theory of the periodic Anderson model in three dimensions, taking 
a cluster of two sites as a basic reference frame. The mean field theory displays the basic features of the Doniach 
phase diagram: a paramagnetic Fermi liquid state, an antiferromagnetic state and a transition between them. 

In contrast with spin density wave theories, the transition is accompanied by a large increase of the effective 
mass everywhere on the Fermi surface and a substantial change of the Fermi surface shape across the transition. 
To understand the nature and the origin of the phases near the transition, we investigate the paramagnetic solution 
underlying the antiferromagnetic state, and identify the transition as a point where the / electrons decouple from 
the conduction electrons undergoing an orbitally selective Mott transition. This point turns out to be intimately 
related to the two impurity Kondo model quantum critical point. In this regime, non local correlations become 
important and result in significant changes in the photoemission spectra and the de Haas-van Alphen frequencies. 
The transition involves considerable / spectral weight transfer from the Fermi level to its immediate vicinity, 
rather than to the Hubbard bands as in single site DMFT 
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I. INTRODUCTION 

The heavy fermion materials have been one of the most 
challenging problems for the condensed matter community in 
the last three decades. The physics of these systems is dom- 
inated by the competition of many interactions, manifesting 
in a variety of phases 1 . Close to the zero-temperature tran- 
sition between two such phases the materials typically show 
non-Fermi liquid behavior on a large temperature range. This 
behavior has been often attributed to an underlying quantum 
critical point (QCP). 

In a subclass of these materials (including CeC\XQ- x Au x , 
YbRh 2 Si 2 and CeTIn 5 , with T=Ir, Co, Rh) it is possible to 
drive the system from an antiferromagnetic (AF) phase to a 
paramagnetic (PM) phase by suppressing continuously the 
Neel temperature via application of pressure, magnetic field 
or doping. The conventional spin-density wave (SDW) the- 
ory of a QCP 2 ,3 ' 4 disagrees with the experimental evidence in 
many of these material a 1 ' 5 ' 6,7,8 . This has triggered a renewed 
interest in microscopic descriptions of the heavy fermion state 
and its possible instabilities to ordered states. More gener- 
ally, not only the critical point but also the formation of the 
heavy fermion state in the regime where the RKKY interac- 
tions are sizeable requires nontrivial techniques for its micro- 
scopic description 9 ^. Heavy fermion materials can be mod- 
eled by the periodic Anderson model which describes a wide 
conduction band of spd electrons and a narrow band of / elec- 
trons which is strongly correlated and hybridized with the spd 
electrons. 

Doniachii established that for large hybridization the / 
electrons combine with the conduction electrons to form a 
heavy Fermi liquid, while for small hybridization the / mo- 
ments order magnetically and hence do not undergo the Kondo 
effect. Thus, changes in the hybridization not only lead to the 
appearance of magnetic long range order but introduce also a 
change in the electronic structure referred to as / — c decou- 
pling or Kondo breakdown (KB). Namely, for large values of 
the hybridization the quasiparticles are composites of / elec- 



trons and conduction c electrons with a Fermi surface which 
obeys Luttinger theorem, counting both the / and c electrons. 
For small hybridization the / electrons decouple from the con- 
duction electron quasiparticles, and the Fermi surface contains 
only the latter. The decoupling between the / and c electrons 
can be viewed as an orbitally selective Mott transition, and 
has been studied with single site DMFT 12 and slave boson 
techniques 13 . As stressed inii, a precise formulation of KB 
requires the elimination of antiferromagnetic long range order 
which breaks the translational symmetry of the original lattice. 
So far this has been achieved by studying models with disor- 
der within DMFT 15 , or by modifying the symmetry group of 
the spin from SU(2) to SU(N)±i. For earlier studies of the KB 
using slave bosons se o 16 i 17 . 

Single site Extended DMFT of various models have been 
carried out. For the Anderson lattice mode l 18 ' 19 , the KB 
and the onset of AF should be viewed as distinct instabilities 
which occur at two different values of the RKKY coupling 
(relative to the Kondo coupling). The same situation occurs in 
the field theory approach of Ref^i. 

On the other hand, extended DMFT studies of a Kondo 
lattice model, in which the conduction electron band is con- 
strained so as to prevent its magnetic polarization, reveal that 
the KB and the onset of magnetism take place at the same 
point in the phase diagra m 20 i 21 ' 22 ' 23 . See however Ref. 2 — . Jus- 
tification for the study of this model is given in Refi^ 5 -. 

In DMFT based approaches the breakdown of Fermi liquid 
theory at the transition is understood in terms of local physics. 
In other works the same breakdown is explained in terms of 
the coupling of the electrons to the long wave length fluctu- 
ations of an emergent gauge field 14 and to the Kondo order 
parameter fluctuation s 1 3,26 . 

In this paper we handle the competition between magnetism 
and formation of a heavy PM state in an unbiased way using 
Cellular Dynamical Mean Field Theory (CDMFT)22. This ap- 
proach improves the treatment of non-local correlations over 
single-site DMFT considering a cluster of sites and provid- 
ing a self-consistency condition on the cluster action. All the 
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correlations whose range is contained within the cluster are 
then treated beyond a simple mean-field level. In this paper 
we consider a cluster of two sites embedded in a three di- 
mensional lattice, which is the minimal unit able to capture 
the competition between PM and AF correlations. A similar 
approach was introduced in an earlier publication 10 based on 
a quantum Monte Carlo impurity solver, but that study was 
limited to large temperatures and small values of U. We over- 
come this limitation by using Exact Diagonalization (ED) as 
an impurity solver 28 . This solver gives access to the ground 
state properties of the system with a finite energy resolution 
(typically 0.005 of the bandwidth in our case). 

In the theory of the Mott transition, the investigation of the 
PM insulating phase has been very fruitful and allowed an un- 
derstanding of the finite temperature state, above the magnetic 
order in many compounds. We pursue a similar approach to 
the periodic Anderson model by investigating both the mag- 
netic and a "parent compound" PM state, to understand better 
the character of the electronic structure. This allows us to de- 
scribe the magnetic state as a result of an additional instability 
of the PM state. 

We show that a phase transition separates the heavy fermion 
paramagnetic phase from an antiferromagnetic phase. Ap- 
proaching the transition from the paramagnetic side the Fermi 
liquid is strongly renormalized. An analysis of the mean field 
PM solution into the AF phase reveals an orbitally selective 
Mott decoupling and localization of the /-electrons. Within 
the resolution of our method, the antiferromagnetic transition 
and the decoupling of the / and c electrons take place for the 
same value of the hybridization. Through a comparison with 
a Wilson's Numerical Renormalization Group (NRG) calcula- 
tion, we further relate the nature of the transition to the prox- 
imity of the system in the parameter space to the quantum 
critical point of the two-impurity Kondo model (2IKM)22i2ii. 

The transition has dramatic consequences on the quasiparti- 
cle spectra and on the topology of the Fermi surface which can 
be detected by photoemission experiments and by Hall effect 
and de Haas-van Alphen frequencies measures. 

The paper is organized as follows. In Sec. |II]we introduce 
the model Hamiltonian and the CDMFT method employed. 
In Sec. [Ill] we analyze the results of the CDMFT calculations 
and compare them with single site DMFT results and Numeri- 
cal Renormalization Group results on the underlying impurity 
problem. Conclusion are presented in Sec. IIVI 



duction electron bandwidth W = I2\t\ = 2 sets the energy 
unit. 

In the large U/V limit the model maps onto the Kondo 
lattice model, where each / orbital is replaced by a spin i 
which is Kondo-coupled to the local conduction electrons with 
a coupling J K = (^p^ + ]u+wp^T\) v2 ' where P is the 
conduction electron density of states at the chemical poten- 
tial. It is also possible to define two energy scales in the 
problem. The first is the single-impurity Kondo temperature 
To = exp(— l/2pJx). It is not known a priori how the "lat- 
tice" Kondo energy is related to To. The second scale is set 
by the magnetic interactions. These are not contained explic- 
itly in this Hamiltonian, but an effective RKKY coupling is 
generated in fourth order perturbation theory in V (second or- 
der in the Kondo coupling, Jrkky ~ J%/W). The RKKY 
coupling oscillates in real space with a period that is related 
to the bare conduction band Fermi momentum, ranging from 
antiferromagnetic at half filling to ferromagnetic in the limit 
of an empty band. A comparison of these two energy scales 
allows to understand the qualitative features of the phase di- 
agram of the modelii. The ratio between Kondo temperature 
and RKKY interaction is an increasing function of V, hence 
for small V we will find a phase in which the magnetic fluc- 
tuations dominate giving rise to an ordered state, while for 
large V the Kondo effect dominates and the / moments are 
screened by the conduction electrons. 

At half filling the Kondo screening of the moments by the 
conduction electrons results in an insulator with a gap induced 
by the hybridization. In this paper we are interested in metallic 
heavy fermions, and we restrict ourselves to study the system 
away from half filling. 

In the path integral formalism the conduction band can be 
integrated out in the action and replaced by an effective long- 
range retarded hopping for / electrons. After these manipula- 
tions the Lagrangian of the system is 



(2) 



where the effect of the conduction band is absorbed in the 
Weiss field 



II. METHODS 
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The three-dimensional periodic Anderson model Hamilto- 
nian is 



H = E c U-a - m E 4^ + v e (/, 



(i) 



and e c 



. cos(g Q ). 



In CDMFT the original three-dimensional lattice is now 
tiled with a superlattice of clusters and an effective Ander- 
son impurity action is derived for a single cluster. Finally a 
self-consistency condition relates the cluster Green's function 
to the local Green's function of the superlattice. The self- 
consistency is more easily written as a matrix form 



where c\ a creates a conduction band electron at site i with spin 
a and f\ does the same for the localized orbitals. The con- 
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where spin index has been suppressed to simplify the notation. 
It is very important to notice that in the large U/V limit the 
impurity action becomes equivalent to the action of a two im- 
purity Kondo model away from particle-hole (p-h) symmetry. 
At p-h symmetry this model is known to have a QCP separat- 
ing a Kondo screened phase from a phase in which the impuri- 
ties are decoupled from the conduction band 29 ' 30 . Away from 
p-h symmetry the QCP is not accessible but still influences the 
physics in an intermediate energy rang o 31 ' 32 . 

In order to solve the cluster-impurity problem we employ 
ED. In this method the cluster action is recast in the form of an 
Anderson impurity Hamiltonian where the continuous degrees 
of freedom of the free electron bath are parametrized by a 
finite number Nf, of sites: 



Hed = £ fl (E f - fi) f ia + U ^2 4™U < 5 > 

,2 i=l,2 



i=l,2 



l=l,N b 



i,l 



In this representation the hybridization-function A M „ = 
— (iuj + fj, — Ef) 6^, related to the non interacting im- 
purity Green's function Q^ v = —{fufv)u—o of the cluster 
impurity problem, becomes 



E 

1=1, N h 



V* la Vua 



(6) 



where the indices i, j = 1,2 represent the cluster sites. In the 
figures presented in this work we used 8 sites in the bath, but 
the results have been tested in few points using 10 sites. 

The ground state and Green's functions of this Hamilto- 
nian are determined via the Lanczos procedure and the self- 
consistency equation in turns allows to derive a new set of 
bath parameters. The self-consistency condition is used to 
determine the A function on a finite mesh of points u> n — 
(2n + 1)tt/ (3. This implies a finite resolution in frequency. 
We used (3 = 100. The process is iterated until convergence 
is reached. To study the interplay of antiferromagnetism and 
Kondo screening we restrict the study to the metallic regime 
close to half-filling where the RKKY interaction is predomi- 
nantly antiferromagnetic. This is achieved by fixing the values 
U = 10, Ef = —5.5 and [i = 0.2 throughout the paper. Being 
a mean field theory, DMFT gives the possibility to selectively 
allow some instabilities and forbid others. Thus we will allow 
the development of AF order and study the realistic material 
phase diagram. However we will also force the system to re- 
main in a PM state to study the "normal state" underlying the 
AF phase. This will prove important to study the KB point 
and to understand the origin of the transitions in the system. 
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FIG. 1: (Color online) Magnetization M and average of the quasipar- 
ticle residue Z eigenvalues for the AF (both for majority and minor- 
ity spin) and PM solutions as a function of V. Inset: Ratio between 
lattice and single-impurity Kondo energies in DMFT and CDMFT. 



HI. RESULTS 

We compute the staggered magnetization M = (n^) — 

(n{|) = — ({n^) — { n 2[)) as a function of V (FigJTJ. For 
small V the RKKY interaction is dominating and the solution 
describes an AF state with a magnetization close to one. In the 
opposite limit of large V, where the Kondo screening is more 
effective, the system is a PM metal. Increasing V in the AF 
phase, the magnetization decreases smoothly and approaches 
zero at a value V* ~ 0.585 where a phase transition occurs. 
Above this critical value the system becomes paramagnetic. 
The value of M at the transition jumps from ~ 0.22 for V = 
0.58 to ~ 0.006 for V = 0.59 and in between the convergence 
becomes extremely slow. The numerical accuracy does not 
allow us to exclude a weakly first order transition scenario, 
although in this case the coexistence region would be very 
small. 

To clarify the nature of the transition we investigate the 
cluster quasiparticle residue Z = (1 — d'E/dui)~ 1 \ L j = o (in 
this case, a matrix in spin and cluster indices). The AF three- 
dimensional SDW theory predicts that antiferromagnetism de- 
velops from a well defined Fermi liquid state and hence a finite 
quasiparticle residue Z at the transition. On the contrary KB 
scenarios predict that Z goes to zero at the transition. In Fig. 
[T]we plot the Z eigenvalues. In the AF phase the two curves 
correspond to minority and majority species. In the PM phase 
we plot only the average of the eigenvalues since the differ- 
ence is not appreciable on the scale of the plot. In both phases 
the eigenvalues of the Z matrix decrease as we approach the 
transition. Precisely at the transition one eigenvalue goes to 
zero, suggesting that the KB is indeed the correct scenario in 
this model. 

Approaching the transition from the PM phase by reducing 
V, we observe a strong reduction of the Z eigenvalues leading 
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FIG. 2: (Color online) Spectral function of the / electrons on the 
two sides of the transition with AF order allowed (upper panel) and 
forbidden (lower panel). The two Hubbard bands are outside the 
energy range. 



to a large effective mass enhancement. This result suggests 
that increasing the correlations by decreasing V from the PM 
side, the Fermi liquid becomes heavier and heavier, and even- 
tually it is more favorable to get rid of the entropy by develop- 
ing AF long range order. The fundamental difference with the 
SDW picture is that here the mechanism responsible for the 
instability is short-ranged in space. The cluster self-energies 
are in fact approximately local on the PM side of the transi- 
tion. This may lead to think that the results are essentially 
similar to a single site DMFT solution of the periodic Ander- 
son model. However there are essential differences. Indeed, 
even setting the off-diagonal self-energy £12 = the CDMFT 
self-consistency equations do not reduce to the corresponding 
single-site DMFT equations. To demonstrate this we com- 
pare in the inset of FigQ]the ratio between the lattice Kondo 
temperature obtained through the relation Tk = A(0)Z and 
the corresponding single-impurity energy To in the two cases. 
While the DMFT result gives a smooth and finite Tk across 
V*, in CDMFT Tk appears to approach a zero value at the 
transition. 

We take advantage now of the possibility of suppressing 
the AF long range order to study the normal PM solution in 
the AF phase. For V < V* this solution is obviously unsta- 
ble to antiferromagnetism and it corresponds physically to a 
highly frustrated system. The analysis of the spectral function 
along the two solutions (allowing and forbidding magnetic or- 
dering, FigEJ reveals that the transition takes place also in 
absence of antiferromagnetism. Namely, at small V, a phase 
with no spectral weight at the Fermi level is realized in both 
cases. The Kondo peak present in the PM phase disappears 
in going across the transition, leaving behind a pseudogap, 
and the / electrons effectively decouple from the conduction 
band, at least within the energy resolution of our calculations. 
The transition is an orbitally selective Mott transition for the 
/ band. This implies that the / electrons do not contribute 
to the Fermi surface for V < V* . This is in marked con- 




FIG. 3: (Color online) Imaginary part of the even self-energy eigen- 
value on the imaginary axis for different values of the hybridization. 
Inset: odd eigenvalue (in a different range of hybridization values). 



trast with the paramagnetic single site DMFT solution where 
there is only a Fermi liquid Kondo screened phase and no sign 
of such a transition or of a pseudogapped phase. Notice also 
that in the CDMFT treatment the / transfer of spectral weight 
takes place from the Fermi level to its immediate vicinity as 
the hybridization is reduced. This should be contrasted with 
the single site DMFT description of the Mott transition in the 
Hubbard model where the spectral weight is transferred from 
the Fermi level directly to the Hubbard bands at much higher 
energies. 

To clarify the Mott character of the transition we plot 
in Fig [3] the imaginary frequency behavior of the two E- 
eigenvalues within the PM solution, named respectively 
E et)e „ = En + E12 and E otW = En - E i2 . Approaching the 
transition from the PM phase the slope at zero frequency of 
both eigenvalues increases, signaling the growth of the effec- 
tive mass in the Fermi liquid. At the transition the even eigen- 
value develops a divergence signaling the end of the Fermi 
liquid description and the beginning of the Mott insulating 
phase. The odd eigenvalue has a similar divergence at a lower 
value of the hybridization as shown in the inset. A qualitative 
understanding of this behavior can be proposed considering 
that E e „ e „ and £ dd represent in a loose sense the lattice self- 
energy calculated around momenta (0, 0, 0) and (tt, it, it) re- 
spectively. Hence the Mott transition corresponds to the first 
appearance of a surface poles in the self-energy around the 
origin of the Brillouin zone. Decreasing further the hybridiza- 
tion the surface of poles moves towards the corner of the zone 
and finally disappears for V ~ 0.43. 

While for V > V * the self-energy was approximately local, 
in the PM solution the development of magnetism for V < V * 
is replaced by a dramatic increase of non-local correlations as 
shown in Fig|4l where we plot the real part of E12 at the lowest 
available frequency in the imaginary axis. 

The behavior of Tk, of the spectral function and of the 
off-diagonal self-energy near the transition is very similar 
to the corresponding behavior found in the 2IKM close to 
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FIG. 4: (Color online) Non local correlations (ReEia, calculated at FIG. 5: (Color online) / and c electron occupation per site. The 
the lowest available imaginary frequency) along the AF and PM solu- transition is marked by the vertical dashed line, 
tions. The vertical dashed line marks the transition. For comparison 
with the PM case the NRG results (without self-consistency). 



its QC P 31,32 ' 33 . To establish this connection we used Wil- 
son's Numerical Renormalization Group to solve the impurity 
model underlying the self-consistency in the PM case for the 
same set of parameters. The difference with the ED calcula- 
tion is that we used a constant density of states for the NRG 
bath, so the results are non self-consistent. The results are nev- 
ertheless qualitatively similar. The main difference between 
the results of the two calculations is that the transition found 
in CDMFT is replaced in NRG by a crossover. We ascribe this 
effect to the lack of a self-consistent bath in the NRG calcula- 
tion. In Fig|4]we compare the non local correlations, finding a 
very good agreement. Other features like the pseudogap for- 
mation around V * are also reproduced. The similarity of these 
two results suggests that most of the physics relevant to the lat- 
tice is already captured at the level of the two impurity prob- 
lem. In this case the onset of magnetism in the lattice model 
can be understood as arising from the vicinity to the impurity 
QCP. Indeed at the 2IKM critical point there are three rele- 
vant operators that correspond to possible symmetry breaking 
instabilities. Susceptibilities in the antiferromagnetic, super- 
conducting and p-h symmetry breaking channels are diverging 
at the impurity QCP. Since p-h symmetry is already broken in 
the lattice system, only the other two instabilities are left to 
provide a way to quench the residual entropy associated to the 
QCP. In this work we did not consider the superconducting in- 
stability but we expect it to play a relevant role (and possibly 
compete with antiferromagnetism) depending on the specific 
details of the band structure of the realistic material. 

Notice however several important differences with the pic- 
ture of local quantum criticality of Ref.— . The reference sys- 
tem describing the physics of the lattice model requires a clus- 
ter of two sites, rather than one, embedded in a medium. The 
critical point of this two impurity model, has two relevant di- 
rections corresponding to varying the hybridization and break- 
ing p-h symmetry. So the impurity model without the self- 
consistency condition has a multicritical point rather than a 



critical point. Generic problems, as the one considered in this 
paper, are p-h asymmetric, and therefore in these cases the 
two impurity QCP is not reachable. However the analysis of 
the NRG low-energy spectrum confirms that in going from 
one side to the other of the crossover (correspondingly, from 
one side to the other of the transition in the CDMFT case) 
the system evolves from a state in which the impurities are 
screened by the conduction electrons (for large V) to a state 
in which the impurities take advantage of the RKKY interac- 
tion to form a singlet and decouple from the conduction band 
(for small V). Hence the low-energy picture for V < V* is 
that of a conduction band asymptotically decoupled from a 
collection of spins arranged in singlets. In the lattice model 
the self-consistency condition transforms the crossover of the 
impurity model into a true phase transition. 

Following the paramagnetic solution we can also observe a 
violation of the Luttinger theorem for V < V* in the form 
of a mismatch between the number of particles and the vol- 
ume of the Fermi surface. This is clear from the plot of the 
number of c and / electrons per site as a function of V (Fig. 
[5). In the PM phase both the conduction and the / bands 
are doped (nf,n c ^ 1). In this case the Fermi surface has 
a volume equal to the doping 6 = n c + n/ — 2. The total 
occupation is not constant because we fixed the chemical po- 
tential rather than the density. For V < V* the occupation 
in both bands becomes essentially constant, implying that V 
is no more effective. In this regime the only contribution to 
the Fermi surface comes from the c electrons which, being 
non-interacting and decoupled from the / band, have a Fermi 
surface with volume n c = 1 + S. These results once again 
support the orbitally selective Mott transition picture of the 
KB 12,13 . Namely, an /-like fluid that becomes incompressible 
beyond the KB point. 

In Fig|6] we plot the fc-resolved total spectral density 
— ilm(G c (k, u))+G{ (k, oj)) at low energies along the param- 
agnetic solution in the two phases. We intentionally consider 
points far from the transition to show that the results are robust 
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FIG. 6: (Color online) Total spectral function as a function of mo- 
mentum and energy on the path connecting the origin T = (0, 0, 0) 
and the corner n = (n, n, n) of the Brillouin zone, (a) V = 0.65 
PM phase, (b) V = 0.30 PM solution in the AF phase. 
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FIG. 7: (Color online) Total spectral function as a function of mo- 
mentum and energy along the same path of Fig. [6]for the AF solution. 
(V = 0.50) 



and do not depend on the vicinity to special points in the phase 
diagram. To obtain this quantity we reconstruct the rotational 
and translational lattice symmetry, broken by the cluster-tiling 
procedure, by periodization of the cumulants 34 . This scheme 
is more appropriate when correlations are stronger. Indeed the 
simple periodization of the self-energy does not reproduce the 
local quantities correctly upon integration over momenta, in- 
troducing spurious states in gapped regions of the spectrum. 
In panel (a) the system is in the PM phase {V — 0.65). It 
is possible to track the origin of the small value of Z in the 
narrow band crossing the chemical potential. The Fermi sur- 
face encloses a "small" region of the Brillouin zone centered 
around the origin. In panel (b) we show the PM solution in 
the AF phase (V — 0.30). It is clear that the volume of the 
Fermi surface in this phase has increased. The difference be- 
tween the Fermi surface volume and the density is equal to 
one within the numerical precision, confirming the violation 
of the Luttinger theorem as expected in the orbital selective 
Mott transition 13 . It is also possible to track the non-dispersive 
/-band at around to ~ —0.2. 

Since the PM solution for V < V* is unstable to antifer- 
romagnetism, to make contacts with experiments we plot in 
Fig|7]the spectral density for the AF solution (V = 0.50). The 
bands are doubled because of the AF ordering. These results 
are representative for the whole AF phase down to small value 
of V. The striking difference with respect to the PM phase is 
the big rearrangement of the spectral weight: now the occu- 
pied region is centered around the ( -| , ^ , -| ) point. This result 
implies a sudden jump across the antiferromagnetic transition 
in all the properties connected with the topology of the Fermi 
surface and it is consistent with the experimental observation 
of a discontinuous behavior of the Hall coefficient 35 and of the 
de Haas-van Alphen frequencies 36 . The contribution of the / 
electrons in this phase is manifest in the two non-dispersing 
bands at uj ~ —0.1 and to ~ 0.45, corresponding to the two 



peaks in Figj2] while the bands crossing the chemical potential 
have mostly c character. 



IV. CONCLUSIONS 

To summarize, we have constructed and investigated a clus- 
ter dynamical mean field theory for the periodic Anderson 
model using a link as a reference frame to capture the competi- 
tion of the Kondo effect and the RKKY interaction. The mean 
field phase diagram away from half-filling features a transition 
separating a paramagnetic heavy electron phase from an anti- 
ferromagnetic phase. In the PM phase, the phase transition 
is accompanied by a strong renormalization of quasiparticles, 
marked by a large enhancement of the effective mass on all 
the Fermi surface. 

By constraining a mean-field metastable paramagnetic so- 
lution into the AF region of the phase diagram, we have been 
able to show that the the transition can be described as an 
orbitally selective Mott transition. The decoupling of the /- 
band from the conduction band marks a localization of the 
/-electrons, as shown by the opening of a pseudogap in the 
single particle spectrum and the appearance of poles in the 
self-energy at the Fermi level, a clear fingerprint of "Mot- 
tness". This is further supported by the incompressibility of 
the system, which presents a constant density by varying the 
hybridization, and by the violation of the Luttinger theorem 
for V < V* along the PM solution. It is important to compare 
this transition, which presents a orbital-selective Mott charac- 
ter, with the traditional DMFT Mott transition of the single 
band Hubbard model. In the present case the weight is shifted 
from the chemical potential to the immediate vicinity form- 
ing a pseudogap rather than being transferred to the Hubbard 
bands at much higher energy. This feature is the product of 
our link-DMFT analysis in three dimensions, strikingly dif- 
ferent from the infinite dimensional DMFT result where the 
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physics is purely local. It would be important in future DMFT 
investigations to better understand the nature of this Mott-like 
transition by studying its cluster-size dependence. 

In the SDW scenario, the fermions which are not coupled 
by the ordering vector are unaffected by the transition. On the 
contrary in our CDMFT study, the transition is accompanied 
by a dramatic rearrangement of the Fermi surface, in agree- 
ment with recent experiment s 35 ' 36 . 

Our findings can be understood in terms of proximity to a 
decoupling QCP with the same properties of the QCP found 
in the 2IKM, namely the impurity model underlying the self- 
consistency equations. Although the lattice system cannot 
reach that critical point due to the particle-hole symmetry 
breaking, the region in which the effect of the critical point 
is relevant is large enough to trigger the lattice antiferromag- 
netic instability 3 1 . If we follow the paramagnetic solution the 
vicinity to the critical point results instead in an increase of 
the non-local correlations £12 in agreement with the behavior 
of the 2IKM. In real systems this effect is preempted by the 
lattice antiferromagnetic instability. It would be interesting 
to realize heavy fermion systems in highly frustrated lattices 
to suppress the antiferromagnetic instability. This would pro- 
vide a physical realization of the paramagnetic solution of the 
CDMFT equations. 

In this framework the antiferromagnetic transition is the 



result of the response of the system to the instability of the 
nearby decoupling QCP, rather than being an intrinsic fea- 
ture. Notice that other correlations are diverging at the 2IKM 
critical point. We expect that depending on the details of the 
model other instabilities (mainly superconductivity) compete 
with antiferromagnetism to order the system and quench the 
residual entropy of the unstable QCP. 

We believe that our description captures the correct physics 
of the periodic Anderson model in an intermediate asymptotic 
regime of energies and proximity to the transition. At higher 
energies, the single site DMFT description is regained. At 
very low energies and very close to the transition, further non 
local effects could turn out to be very important. We do not 
know if our treatment fits within the local quantum criticality 
paradigm where no spatial anomalous dimensions are gener- 
ated. To elucidate these questions would require the investiga- 
tion of the dynamical susceptibility beyond the cluster DMFT 
approximation, which is beyond the scope of this work. 

We would like to thank M. Fabrizio, A. Georges, F. 
Steglich, P. Coleman, C. Pepin, I. Paul, O. Parcollet, P. Sun, 
K. Haule, E. Lebanon, M. Capone, T. Stanescu, J.C. Domenge 
and L. de' Medici for helpful discussions and comments. This 
work was supported by the NSF under Grant No. DMR 
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